LEAD  I  DISCOVER  I  DEVELOP  I  DELIVER  ^ 

Theoretical  studies  of  nanociuster 

formation 

26  May  2016  _ 


Boatz 

Principal  Research  Chemist 
Aerospace  Systems  Directorate,  RQRP 
Air  Force  Research  Laboratory 

This  briefing  contains  information  up  to: 


Distribution  A:  Approved  for  public  reiease;  Distribution  Unlimited.  PA  Clearance  No:  16215 
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1.  Introduction 

-  background,  technical  approach 


2.  Core-shell  nanoclusters  (Mg/Cu,  Si/AI,  etc.) 

-  energetic  additives  for  propellants,  explosives 

-  gas  generators 

-  biocidal  defeat  agents 

3.  Summary  and  conclusions 
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Core-shell  nanoclusters  such  as  SiAl^,  Ni^AI^,  AI„(CuO)^,  etc.  may  be  useful 
ingredients  in  propeilants  and  expiosives 

-  higher  energy  densities  than  organics  (~  3x  RDX) 

_ -  some  are  resistant  to  surface  oxidation  (i.e.,  “magic  clusters”) 


Helium  droplet 
experiments  at 
AFRL/RW 


Can  core-shell  nanoclusters  be  formed  under  cryogenic  conditions  (i.e.,  in 
helium  dropiet  experiments)  via  stepwise  condensation;  i.e.,  what  are  the 
energy  barriers  (if  any)  to  stepwise  addition  of  atomic  Ai? 

_ SiAI„  +  Ai  ^  SiAI„^, ,  [Aln]-  +  Ai  ^  [AI]- _ 
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Cu^Mgy  core-shell  nanocluster 

inversion 


In  helium  droplet  experiments,  Mg  atoms  were  captured  in  first  pickup 
cell,  followed  by  capture  of  Cu  atoms  to  form  Cu^Mgy  core-shell 
nanoclusters.  However,  scanning  transmission  electron  microscopy 
(STEM)  measures  show  cluster  inversion  occurred  to  produce  MgyCu^C!) 


a)  copper  atoms 

b)  magnesium  atoms 

c)  oxygen  atoms 

d)  composite  image 


MP2/ACCT 


R=  3.780  A 
E.=  2.65kcal  mol*i 


R=  2.396  A 
E.=  19.73kcalmo|-^ 


R=  2.236  A 

B.E.=  46.39  kcal  moK 


HCAnfstrom) 
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1 .  structure  of  Mggo  cluster 
was  fully  optimized. 

2.  Two  Cu  atoms  were 
placed  on  opposite  sides  of 
Mggo  and  structure 
reoptimized. 

3.  Distance  between  Cu 
atoms  was  decreased  in 
steps  of  0.25  A,  held  fixed, 
and  remaining  DOF 
reoptimized. 

4.  Total  energy  plotted  as 
function  of  fixed  Cu-Cu 
distance. 
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r  (Cu-Cu)  (A) 
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5 _ 9 _ 10 _ 

Energy  barrier  is 
~0.5  kcal/mol 


DFT  calculations:  B3PW91/aug-cc-pa)CVTZ(-PP)  level 
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Drop  in  energy  in  Cu2Mg3o  is  not  due  to  Mg3o  rearranging  to  more  stable  structure. 
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Cu3oMg2 
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Energy  barrier  is 
~6  kcal/mol 
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Drop  in  energy  in  Mg2Cu3o  is  not  due  to  CU30  rearranging  to  more  stable  structure. 


50 


-40 


-■-Cu30 


3  4  5  6  7  8  9 


r  (Mg-Mg)  (A) 
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Formation  of  SiAl 


MRMP2(10e,12o)/aug-cc-pvtz  level  of  theory 

Cyclic  isomers 

5Ai,  -0.63  eV 


3A”  ,  -1 .80  eV 


iAi ,  -2.00  eV  3b^  _  .2.04  eV 
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Formation  of  S/A/2 


MRMP2(10e,12o)/aug-cc-pvtz  level  of  theory 


Linear  isomers 


in„.-1-28eV 
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Summary  and  Conclusions 


•  Mg/Cu  core-shell  nanoclusters 

•  Helium  droplet  experiments  show  inversion  of  Cu^Mgy  clusters  to  MgyCu^. 

•  Cu  atoms  diffusing  into  Mg3o  and  vice-versa,  have  been  modeled  using  DFT. 

•  Estimated  barrier  for  Cu  atoms  to  migrate  into  Mg^  is  <  1  kcal/mol. 

•  Estimated  barrier  for  Mg  atoms  to  migrate  into  CUp  is  6  kcal/mol. 

•  Calculations  are  consistent  with  observed  Cu/Mg  inversion. 


•  SiAl^  clusters 

•  SiAl2  has  multiple  local  minima  which  are  more  stable  than  ground  state  SiAl 
m  +  Al  (2P) 

•  Cyclic  (C2V  and  Cg)  and  linear  ^^^l  C^ov) 

•  Singlet,  triplet,  quintet  states 

•  At  long  SiAl - Al  separations,  preferred  approach  is  linear 

•  AI-Si-AI:  ^riu  without  a  barrier 

•  Si-AI-AI:  can  form  without  a  barrier 

•  Barriers  for  linear  ^  cyclic  isomerizations  TBD 

•  In  helium  droplet  environment,  linear  quintet  states  may  be  formed. 
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Mg^  benchmark  calculations 


Mgp  clusters 

“Closed  shell”  atomic  configuration  [(1s)2(2s)2(2p)®(3s)2]  suggests  that 
weak  dispersion  interactions  will  be  important.  Need  to  consider 

-  core-core  and  core-valence  correlation 

-  correlation  method  (MP2,  CC,  DFT) 

-  “active”  electrons  to  be  correlated  in  MP2,  CC 

-  suitable  DFT  functional  for  larger  Mg^ 

clusters  (up  to  n  «  100) 


Method 

cc-pwCVDZ 

cc-pwCVTZ 

cc-pwCVQZ 

MP2 

23.1  /  3.042 

28.5  /  3.013 

29.1  /  3.011 

CCSD(T) 

16.7  /  3.100 

tbd  /  3.064 

tbd  /  3.065 

DFT/B3PW91 

26.5  /  3.092 

26.5  /  3.091 

26.4  /  3.092 

DFT/PBE 

34.6  /  3.070 

34.5  /  3.070 

34.3  /  3.070 

DFT/PBEO 

31.5  /  3.078 

31.5  /  3.078 

31.4  /  3.078 

DFT/M06 

30.8  /  3.028 

30.4  /  3.025 

TBD  /  TBD 

DFT/M11 

19.2  /  3.134 

TBD  /  TBD 

TBD  /  TBD 

Calculated  binding  energies  used  to  determine  size  of  helium  droplet  needed  for  evaporative  cooling 
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Calculate  potential  energy  surfaces  of  stepwise  atomic  addition  reactions; 
_ e.g.,  SiAi„  +  Al  ^  SiAI„+i _ 


Diatomic:  Si  pP)  +  Al  pP)  ^  SiAl  ) 


Triatomic:  SiAl  )  +  Ai  (^P)  ^  SiAig  (triplet  or  quintet) 


302 ,  -2.67  eV 
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/■ 


SM/  +  4/  ^  SM/. 


SiAl  +  Al  (2P) 


R(AI-AI) 


AlSi  (^X)  +  Al  (2P) 


R(Si-AI) 


^ 


m) 


^ ^ 


(^n) 


MRMP2(4e,5o)/aug-cc-pVTZ 
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